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Abstract 



The interaction of quarkonium with nuclei is studied in the wiq — > oo limit of QCD, 



where the binding energy is found to be exactly computable. The dominant contribution 
to the interaction is from two-gluon operators. The forward matrix elements of these two- 



gluon operators can be determined from the QCD scale anomaly, and from deep inelastic 
scattering. We apply our results to the T and J/i(j, treating the QQ interaction as purely 
Coulombic. We find the T binds in nuclear matter with a binding energy of a few MeV, 
while for the J/ip binding is of order 10 MeV. For the J/ijj in particular we expect 
confinement effects to produce large corrections to this result. 
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There has been some recent interest in the possibility of bound states occurring be- 
tween quarkonium and nuclei |l] 0- By scaling -KT- nucleon scattering from high energies 
and modifying the interaction due to the reduced size of charmonium the binding energy 
of the r\ c in nuclear matter [|J is estimated to be ~ 30 MeV. The motivation for study- 
ing the interactions between quarkonium and nuclei is clear. The interactions of nuclei 
cannot be computed starting directly from the QCD Lagrangian. However, they can be 
parameterised by using a phenomenological Lagrangian with nucleon and meson degrees 
of freedom. Thus a study of nuclear forces provides some understanding of the dynamics 
of nucleons and the pions which are produced by the spontaneous breaking of chiral sym- 
metry. The interactions of quarkonium allows us to study a different aspect of QCD, the 
interactions of gluons. The most significant difference between the nucleon- nucleon poten- 
tial and a quarkonium-nucleon potential is that there is no valence light quark exchange 
in the latter. The coupling constants of the interactions between quarkonium with pions, 
nucleons and other light states arise solely from gluon induced interactions. 

In this paper, we study the interactions of a QQ bound state with nuclear mat- 
ter, in the limit that the mass uiq of the quark, and therefore the inverse radius 
Tq 1 ~ a s yoj m Q °f the QQ bound state, is much larger than the QCD scale Aqcd- 
We then apply this to the cases of experimental interest, the cc and bb bound states, and 
calculate their binding energies in nuclear matter. 

At scales \x <C r^ 1 , the QQ bound state looks like a colour singlet, pointlike object. 
Its dynamics will be described at low momenta by an effective Lagrangian in which non- 
renormalizable terms are suppressed by powers of the compositeness scale Aq = Tq 1 (~ 
a(AQ)rriQ for a Coulomb bound state). The gauge invariant gluon couplings of lowest 
dimension are the dimension 7 operators^ 

= E T3 ( P v p v ~ V tvK) (ceO e + C B B ) , 

O e = -G» aA G» a A v„v u , O b = ^G af3A G A p - G^ A G vA v^v vl 

where P v and create pseudoscalar [r\ c or r]b) and vector (J/ip or T) mesons respectively, 
with four velocity t> M . All the sensitivity to the detailed structure of the meson is encoded 
in the coupling constants eg and eg, which must be calculated from the underlying QCD 
dynamics. 

1 This is a just a covariant formulation of the leading terms in the usual multipole expansion S. 
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Some explanation of (|l|) is in order. Since the interaction of a heavy state with light 
degrees of freedom does not change the velocity of the heavy state in the m,Q — > oo limit, 
the velocity is conserved, and it is convenient to describe mesons with different velocities 
by distinct fields[|J. The fields Vtf and P v are related to the standard fields P and by 
the field redefinitions 

V£(x) = V2me imVaXa V ll (x), P v (x) = V2^e imVaXa P(x), (2) 

where m is the meson mass, and physical vector states satisfy v^V^ = 0. Writing the 
momentum as 

pf* = mv» + , (3) 

we see that derivatives acting on the redefined fields bring down factors of the small "resid- 
ual" momentum /c M , so the derivative expansion does not contain potentially troublesome 
terms proportional to tuv^/Kq. Also, in the heavy quark limit, the quark spin decouples 
and the resulting spin symmetry [|J forces the masses and couplings of the pseudoscalar 
and vector bound states to be equal (the minus sign in ([]]) is because a physical vector 
particle has spacelike polarization). There is also no term which converts a pseudoscalar 
to a vector at leading order in 1/mg. Hence, in the uiq — > oo limit, the binding to nuclei 
of a vector meson is equal to that of the corresponding pseudoscalar meson. 

In the rest frame of the QQ bound state, v = (1,0,0,0), and the operators Oe and 
Ob become 

O e = E A ■ E A , O b = B A ■ B A , (4) 

where E A and B A are the colour electric and magnetic fields. (This explains the choice 
of operators and coefficients in Eq. fll]).) The coefficients ce and cb can be obtained by 
computing the energy shift of the QQ bound state in an external electric and magnetic field, 
i.e. by computing the quadratic Stark and Zeeman energies. The dominant contribution is 
the quadratic Stark effect, since magnetic effects are suppressed by a factor of v/c ~ ck s (Aq) 
relative to electric effects for each replacement E — >• B. 

The operators Oe and Ob can be written as linear combinations of the gluon energy- 
momentum tensor, 

T Zon = \^ v G a0A G% - G^ aA G^ A , (5) 
3 



and the divergence of the scale current,! 

3% = T% = ^G^ A G%. (6) 

Note that the energy-momentum tensor T^ v that occurs in Eq. (Q) is the full energy- 
momentum tensor, whereas Tj^ on in Eq. (|5|) is only the gluon contribution to the energy- 
momentum tensor. The expressions for Oe arid Ob are 

sr\ rp\lV 9 rpOL 

U E — i gluon t, M t '^ ~~ 23(a) ' 

g (7) 



The operators (9^; and Ob in the effective Lagrangian Eq. ([!]) are renormalised at the scale 
Aq, and the operator T^ on in Eq. (0) is also renormalised at Aq. (The operator is 
scale independent.) The QCD beta function at the scale Aq can be written as 



^9 = P(9) = -b Q ^ + ..., (8) 



so that Eq. (|7]) can be rewritten as 



(9) 



Ob = - j^^y 

The gluon operator in Eq. ([]]) can be rewritten in the form 

ceO e + c b Ob = (c E + c B ) (t^v^) + (ce - c B ) ^ A ^ (10) 

The matrix elements of T^ on and T£ in a nucleon state are known at zero momentum 
transfer, 

<p|T Q »=2M 2 , 

(P\ ^gluon |P> = - I^V) , 

where M is the nucleon mass. The value of V2, the gluon momentum fraction in the nucleon, 
is measured in deep inelastic scattering, and its value can be extracted from the F\ (or F 2 ) 



2 For simplicity, we will work in the limit that the light quark masses are zero. It is straight- 
forward to include light quark mass effects Q. The trick of using the scale anomaly to compute 
the matrix element of G 2 is due to Voloshin and Zakharov[Ml. 
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structure function. The usual value quoted in the literature is V^/U 2 = 16 GeV 2 ) = 0.44 
. The fi dependence of V2 may be found in || . 

The forward scattering amplitude of a QQ meson off nucleons is completely determined 
in terms of the measured value of V2, 

Mm _ mAQ) (^) M * (7 . _ J, + 2M , ) (12) 

where 7 (= v ) is the time dilation factor in the rest frame of the nucleon. The forward 
scattering amplitude is proportional to the energy shift of the QQ meson due to its inter- 
action with nuclear matter. The binding energy is given in terms of M. (neglecting Fermi 
motion) by 

Mn 

^binding = ^77 ' \^) 

where n is the number of nucleons per unit volume in nuclear matter. The factor of 2M 
divides out the standard relativistic normalisation of nucleon states used in computing Ai 
in Eq. (|T2|). The binding energy of QQ states in nuclear matter is 

^ = mM) (7 » _ j) + , ) (14) 

where p = Mn is the mass density of nuclear matter. 

To determine the value of the binding energy of a QQ state, we need expressions for 
ce and cb- The coefficient cb is suppressed relative to ce by a factor of ck s (Aq) 2 , and can 
be neglected. The computation of ce is subtle, because gluon interactions cause mixing 
between the singlet and octet QQ channels. In the heavy quark limit, vtiq — > oo, in which 
the QQ system is a Coulomb bound state, a careful derivation of ce has been given by 
Peskin |J and we will quote his result. For the Is state, he finds 



c E 14 7r^ 3 



A*, 27 



r% (15) 



where N c is the number of colours, and tb is the Bohr radius 

_i 2 



r 



B 



-OLstr^mq. (16) 
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Peskin used the 1/N C expansion in obtaining the expression for ce in Eq. (|T5|). This 
approximation was used merely to give an expression which can be written in a simple 
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form. The error in not keeping the 1/N C corrections was estimated to be approximately 
10%, which is smaller than the a s (Ag) and 1/toq corrections that we have neglected. 
Combining Eqs. ( |T5| ) and gives 



14 3 

^binding = Tfj^BP 



V 2 (A Q ) (7 2 -|)+ 2 " 



6q« s (Aq) 



(17) 



Treating the T as a Coulomb bound state, we find from (^Bp r^ 1 ~ 1.1 GeV. Using a 
more realistic potential, Quigg and Rosner |TD[ calculated (r)y ~ 1.0 GeV -1 , correspond- 
ing to a Bohr radius r^ 1 ~ | (r) _1 ~ 1.5 GeV. To the degree to which our results are 
sensitive to the nonperturbative region of the potential, they are not to be trusted, and 
since t/bmding goes like r|, this introduces a large uncertainty into our results. We give an 
idea of this uncertainty by giving Landing for both values of . For an T or ryj,, we find 
the binding energies for the state at rest in nuclear matter of 

Ending ^ 4 MeV (r~ B l = 1.1 GeV, a s = 0.36), 

(18) 

~ 2 MeV (r" 1 = 1.5 GeV, a s = 0.30). 

Here we have taken p = 1.2 x 10" 3 GeV 4 , and b Q = 27/3. 

For the J/ip, the validity of our expansion becomes markedly more suspect: from (|i~6| ) 
we find r^ 1 ~ 640 MeV, while from |10| we find r^ 1 ~ 750 MeV, where we really do not 
trust our perturbative results. However, proceeding as before we find 

^binding ^ 11 MeV (r" 1 = 640 MeV, a s = 0.6), 

(19) 

~ 8 MeV (r^ 1 = 750 MeV, a s = 0.5). 

Of course we do not believe this result is anything more than a rough estimate for the J/if) 
binding energy, as the corrections to ( |i~9"D will be of order unity. 

Note that our computation unambiguously determines the sign of the energy shift of 
QQ mesons at rest in nuclear matter, at least in the uiq — > oo limit. The spin dependent 
corrections are suppressed relative to Eqs. (^8|) and (|T9| ) by an additional factor of a s , so 
that the binding energy of the r/i, and T are approximately equal. The same arguments 
will also apply to the J/ijj and rj c , but the spin symmetry breaking effects will be larger. 

There are two trivial corrections to the above expression Eq. (|T7j), which we have not 
included because they are small compared to the ck s (Aq) and 1/mg terms that we have 
neglected. The first is that there is a binding energy correction to the relation between 
the nucleon mass density, and the mass density of nuclear matter. The second is that 
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the structure functions of bulk nuclear matter differ from those of nucleons, so that A in 
Eq. ( |T?D should properly refer to the gluon distribution measured using the F2 structure 
function for nuclear matter. 

The results obtained in this paper show that the scattering amplitude of a QQ me- 
son with nuclear matter scales as rg, and the cross-section as Vq (up to the logarithmic 
dependence in a s )$ This does not agree with the assumption in the literature that the 
cross-section scales like the area of the meson. We have also shown that the sign of the 
interaction energy is negative, so that the meson binds to nuclear matter. There are no 
pion corrections to our result, Eq. (|T7j). All pion contributions are included in the matrix 
element of the gluon operators in nuclear matter. The only role of the pions in this cal- 
culation is to modify the density and pressure of nuclear matter. We emphasise that our 
results are rigorous predictions of QCD and are completely model independent in the limit 
mQ — > 00. 

The largest source of error in our computation is the use of Coulomb wavefunctions 
to compute ce- Non-perturbative confinement effects can be included by using a more 
realistic potential model to compute the bound state wavefunction. Unfortunately, the 
computation of ce involves intermediate colour octet states, and thus the QQ Hamiltonian 
in the colour octet channel. In perturbation theory (the computation of Peskin |J), one 
replaces the attractive — 4a s /3r attractive potential in the singlet channel by a repulsive 
a s /6r repulsive potential in the octet channel. The corresponding modification when non- 
perturbative effects are included in the singlet channel is more complicated, and is being 
investigated further. 

The binding energy was computed in Eq. ([17]) for the Is ground state. One can 
compute the binding energies of the excited states using the values of ce given in || for 
the excited states. Since the binding energy goes like rg, the excited states of quarkonium 
will be more tightly bound in the tuq — > 00 limit. Naively applying ( |T7j) to the ip'{2s) 
gives 

^Lling = (jt) ^L g ~ 700 MeV, (20) 

3 The r 3 dependence of the two gauge boson-meson scattering amplitude also shows that the 
Rayleigh scattering cross-section is proportional to r 6 u! 4 . The nuclear binding of QQ mesons arises 
from the same operator that is responsible for Rayleigh Scattering, with the electromagnetic field 
strength tensor replaced by the gluon field strength tensor. 
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suggesting that the ifi'(2s) has a lower energy in nuclear matter than the ground state. 
The large radius of the 2s is clearly an indication that the multipole expansion is breaking 
down, so this result cannot be trusted (but it might still be true). 
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